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We report an investigation of electron conduction in oligophenyl based double barrier molecular 
device. We have carried out analytical calculations and numerical simulations on isolated molecules, 
consisting of aromatic tv conjugated system made up of three phenyl rings separated by insulator 
groups —CH2 — , —S1H2 — , —GeH2— and —SnH2 — . We show analytically as well as numerically 
that when the two insulator groups are different an asymmetric electron transport in the presence 
of external electric field is possible, thus resembling a diode which allows one way electron transport. 



INTRODUCTION 



Q \ Recent progress in experimental and theoretical techniques have triggered enormous possibilities of using organic 

molecules as a fundamental complement to silicon in semiconductors [lM3[- This has opened up possibilities of 
fabricating logical circuits at the single molecule level [4]. Such molecular state electronics requires the precise 
knowledge of electron transport in organic molecules The finite size of atoms, molecules and clusters offers an advantage 
in designing tailor-made materials to perform desired functions by exploiting the quantum mechanical phenomena 
of electrons in a confined space. Besides such exciting developments on the design and fabrication of molecular 
devices, experiments have been performed on the electrical properties of single molecules [5|, and highly sophisticated 
simulation tools for charge transport through molecules have been developed [6-9]. Recently, using a single benzene 
molecule, the first molecular transistor was developed [10]. Several molecular device concepts are based on the 
existence of potential barriers within molecules, which prevent the free flow of electrons and thus create a tunable 
system [Ill4l5|. Aviram and Ratner [16| proposed such a device in 1974 based on an organic molecule made of tt— 
donor and tt— acceptor moieties linked by a saturated spacer and sandwiched between two metallic electrodes. This 
concept has recently been verified experimentally [17j . The presence of potential barriers is essential also in devices 
which are closer to traditional electronic electronic components. 

A diode is a two-terminal switch consisting of source and drain that can turn a current on and off depending on 
the direction of flow. In the case of molecular diodes, switching the current 'on or off' requires tuning the molecular 
discrete energy levels. In an experiment in 1999 [18], it was shown that a single molecule of 2' — amino — 4 — 
ethynylphenyl — 4' — ethynylphenyl — 5'nitro — 1 — benzenethiol could behave like a resonant tunneling diode. The 
resonance effect of tunneling the electrons was observed at 2.1V bias with a peak current of In A. The molecule 
consists of three benzene rings connected to each other by the acetylene group. Two H atoms of the central ring 
C^> J are replaced by the NO2 and NH2 group on opposite sides. This asymmetrical configuration enables the molecule 
susceptible to change its configuration under perturbation. 

In the present work we are interested in a quantitative investigation of molecular structures containing aliphatic 
chains acting as barriers inserted between aromatic tt— conjugated rings, focusing on systems made up of three 
phenyl rings separated by insulator groups — CH2 — , —S1H2— ,— GeH2~ and —SnH2—. We show analytically as well 



as numerically that when the two insulator groups are different an asymmetric electron transport in the presence of 
external electric field is possible, thus resembling a diode which allows one way electron transport. To test the barrier 
properties we have performed analytical quantum mechanical calculations and also performed numerical simulations. 
We have based our analysis on the analytical results obtained for the occupation probability of the energy levels, 
the derealization of the Highest Occupied Molecular Orbital (HOMO) and the HOMO-LUMO energy gap (HLG). 
The HLG is a critical parameter for the molecular admittance because it is a measure of the hardness of the electron 
density [12[. The larger the HLG, the more stable the molecule, and therefore the harder it is to rearrange its electron 
density under the presence of an external electron. Thus the molecule presents lesser admittance to the incoming 
electron. 

THE MODEL AND ANALYTICAL RESULTS 

The source of conductivity for a polypheny lene wire is a set of delocalized tt— type molecular orbit als above and 
below the plane of a single benzene ring. In this section, we present a quantum mechanical analysis of electron 
density behaviour of the polypheny lene wire in which insulator groups — CH2 — , —S1H2—, —GeH2— and —SnH2 — 




FIG. 1: Schematic representation of the system we are investigating. The X and Y are any one of the insulating groups 
—CH2-, —SiH2 — , —GeH2— and —SnH2 — . If X ^ Y then an asymmetry is generated in the electron conduction process 
when an electric field is applied. 

are inserted between two aromatic rings and an external electric field is applied. The system considered is shown in 
Figure 1, where X and Y are the insulator groups —CH2-, —S1H2—1 —GeH.2 — and —SnH2— . If X ^ Y then an 
asymmetry is generated in the electron conduction process when an electric field is applied. We now present here a 
simple three state quantum mechanical model which provides an intuitive understanding of the role played by the 
aliphatic fragments in the triphenyl system. 

The Hamiltonian for a single electron in a system with three levels E±, E2 and £^3, coupled through an external 
perturbation, can be written, in second quantized formalism, as 
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where rii,i = 1,2,3 are the occupation numbers of the three levels accessible to the electron, b \ and b{ are the 
creation and annihilation operators for an electron in the i th level, t\ and £ 2 represents the interlevel coupling or the 
tunneling matrix element. The second term takes into account the perturbing potential V through its value Vi in the 
location corresponding to the i th level. If the Hamiltonian is expanded in the following simple basis 
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|3 >, n\ = 0, n 2 = 0, n 3 = 1. 
and written in matrix form, it reads, in atomic units, as 
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We now define the energies e\ = (E2 + qV2) — (E\ + qV\) and e 2 = (£^3 + ^^3) — (E\ + qV\) and obtain the three 
eigenvalues of the above Hamiltonian as the roots of the following cubic equation 



A 3 - A 2 (ei + e 2 ) + A(eie 2 - t\ - t\) + e 2 t\ = 
Now for the simple case e 2 = and e\ = e, we obtain two eigenvalues as 
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FIG. 2: Occupancy of the first level of the analytically solvable model as a function of the perturbation e, for three different 
values of the ration txjti. 
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In Figure 2 we plot the occupancy of the first level \c\ | 2 as a function of e for three values of the ration t\jti- Here 
c\ is the first coefficient of the eigenvector of the ground state . The analytical expression for c\ is 
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The parameters t\ and £2 describes the effective interaction between the orbitals of the three aromatic rings and 
therefore will be smaller when an aliphatic group breaks the 7r-conjugation and spaces the rings. An asymmetry is 
produced when t\ ^ £2 • In Figure 2, when t\/t2 = 4.0, the single electron will find it easier to tunnel from the first 
ring to the middle ring through the insulator X as compared to the path starting from the middle ring to the third 
ring via insulator Y. Since t\ > £2, the internal barrier between the first ring and the middle ring is small and the 
aromatic fragments are closer together and strongly interacting. On the other hand when £1/^2 = 0.1 exactly opposite 
is true and hence the electron density will tend to localize on the first ring as evident from Figure 2 when |ci| 2 — >• 1 
at a smaller value of e compared to the cases t\/t2 = 1.0 and £1/^2 = 4.0. This asymmetry that is introduced when 
fi/t2 can be utilized to design a tunable molecular tunneling diode. The curve in Figure 2 for the highest value of 
ti/t2 exhibits the desired unidirectional (left to right) flow of electron current if initially the electron is localized on 
the left ring. 



NUMERICAL RESULTS 



All calculations were performed using Hyperchem 7.0 software. To obtain good starting geometries, the initial 
configurations were optimized at the Hartree-Fock (HF) level. The density functional theory (DFT) method was used 
on the optimized geometries obtained from HF calculations. We have used the 6 — 31G* basis and hybrid functional 
B3PW91, which combines the Becke exchange (B) functional, and the Per dew- Wang- 91 correlation functional, both 
of which are nonlocal generalized gradient- approximated functionals. Changes in the molecular orbitals are excellent 




FIG. 3: Optimized structure of the TP — CH2— system with the HOMO distribution (top) and the extend of the LUMO orbitals 
(bottom). Symmetrical distribution of HOMO and LUMO is visible. 




FIG. 4: Shape of the HOMO orbitals of the TP 
atomic unit = 5.142 x 10 n V/ra. 



CH2— system in the presence of electric field E — ±0.001 atomic units. 1 



indicators of the molecular electron transport. It is important to understand the electronic molecular transport under 
the action of an external electric field along the axis of the molecule. In this context the HLG energy and the spatial 
extent of the HOMO are good candidates. We need to quantify two effects regarding the transport of electrons in 
molecules to investigate their electrical properties. First, the HOMO effect needs to be quantified because the charge 
transfer from one end of the molecule to the other end becomes harder to achieve when the HOMO is delocalized over 
the entire molecule. The second effect is due to a directly quantifiable property, the HLG. When this gap decreases, 
the molecular admittance increases. First, we analyze the case when X = Y = —CH2 — . Figure 3 shows the extend 
of the HOMO and LUMO orbitals for the case X = Y = —CH2 — . Clearly we notice the symmetrical nature of the 
extend of HOMO and LUMO across the molecule due to the symmetric nature of electron withdrawing and electron 
donating groups. This is the case of t\ = £2- The geometry optimization was performed initially by the HF and then 




FIG. 5: Optimized structure of the TP — S1H2— system with the HOMO distribution (top) and the extend of the LUMO 
orbitals (bottom). Asymmetrical nature of the LUMO and the HOMO is visible. The distance between the two phenyl rings 
via the —CH2— spacer is 2. 669 A and via the —SHI2— spacer is 3. 136 A. This gives rise to the asymmetrical tunneling matrix 
element. 

by the DFT method. Interestingly the two methods yield different HLG energy. Table 1 shows the comparison of the 
electronic energy levels calculated using HF and DFT methods for the molecule triphenyl system with two —CH2 — 
spacers (TP — CH2—). 



Table 1: Comparison of the electronic energy levels of TP — CH2 — calculated using HF and DFT methods. 



Method 



Eromo (eV) 



Elumo (eV) E LU mo+i (eV) 



HLG(eV) 



DFT 
HF 



-4.73 
-9.12 



2.34 
0.42 



2.39 
0.43 



7.07 
9.54 



A significant difference in the HLG is noticed from Table 1. The reason for such a difference in the results obtained 
from HF and DFT method is due to the absence of electron correlations in the HF method. Electron correlation effects 
are expected to play important role in this molecule. Any incoming electron from the metallic electrodes connected 
to the one end of the molecule is expected to propagate over the LUMO + 1 under a suitable bias voltage required to 
transfer one electron through this molecule. We assume in the first approximation that the Fermi level of the metal 
electrode lies in the middle of the HOMO and LUMO energy levels of the molecule. To determine the bias voltage 
required (E v )to transfer one electron through this molecule is estimated as [14] 



??#« = 



HLG 
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(10) 



Here HLG represents the energy gap between HOMO and LUMO energy levels. Elumo and Elumo+i represent 
the energy for the LUMO orbital and LUMO + 1 orbital, respectively. The E v for the molecule TP — CH2— from 
the above Eqn.(10) is estimated to be 3.59eT/. 

Figure 4 shows the spatial orientation diagram for the HOMO energy levels of the molecule TP — S1H2— (triphenyl 
system with one —CH2— and one —SiH2~ spacer) in the presence of electric field E = ±0.001aw. Due to the 
symmetry of the molecule, the HOMO energy levels are symmetrically displaced around the central phenyl ring. This 
is a consequence of the fact that t\ — £2. 

When the insulator spacer is not identical (I / 7), the situation varies radically. Figure 5, displays the double 
barrier molecular RTD for the case when X = —SiH2— and Y = —CH2— giving rise to t\ 7^ £2 and t\ < ti (The 




FIG. 6: Shape of the HOMO orbitals of the TP — S1H2— system in the presence of electric field E — ±0.001 atomic unit. 1 
atomic unit = 5.142 x 10 11 V/ra. The top figure is for E — 0.001 atomic unit and the bottom figure is for E — —0.001 atomic 
unit. 
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FIG. 7: Plot of the energy gap HLG versus electric field. For the molecule TP — CH2 — , 
absence of any MRTD behaviour. On the other hand, for the second molecule TP — S1H2- 
barrier on the positive electric field side. 



the plot is symmetric indicating the 
- , the electron experiences an energy 



distance between the two phenyl rings via the —CH2— spacer is 2.669A and via the —S1H2— spacer is 3. 136 A This 
gives rise to the asymmetrical tunneling matrix element. A larger distance indicates a smaller tunneling matrix 
element). This asymmetric charge distribution gives rise to asymmetric spatial orientation of the HOMO and LUMO 
energy levels. This is exactly what one requires to design a double barrier molecular resonant tunneling diode. Figure 
6 shows the spatial orientation diagram for the HOMO energy levels of the molecule TP — SHI2— in the presence 
of electric field E = ±0.001aw. Due to the asymmetry of the molecule, the HOMO energy levels are asymmetrically 
displaced around the central phenyl ring. This is a consequence of the fact that t\ ^ £2- For E = O.OOlaw the HOMO 
is delocalized over the entire molecule indicating a reduction in the molecular admittance. In this case, the shift of 
the electron density is obviously less pronounced, since all the n- electrons are strongly involved in chemical bonds 
and transfer of an electron from one ring to the other becomes clearly much harder to achieve. On the other hand 
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when E = —O.OOlau the HOMO energy levels are localized on the central and the right phenyl rings which indicates 
localization of the excess electron on the right phenyl ring. The E v for the molecule TP — S1H2— in the absence of 
electric field is estimated to be 3.25eF which is lower than that of the previous molecule TP — CH2 — (3.59ei>). A 
plot of energy gap HLG versus electric field is shown in Fig. 7. For the molecule TP — CH2 — , the plot is symmetric 
indicating the absence of any MRTD behaviour. On the other hand, for the second molecule TP — SiH2~, the electron 
experiences a energy barrier on the positive electric field side. This basically indicate that on application of an electric 
field in the positive direction of the molecular axis, the electron experiences a higher barrier as compared to the case 
of an electric field in the negative direction of the molecular axis. This asymmetric response of the molecule to the 
external field is useful for possible molecular tunneling diode application. Replacing Si with Ge in the TP — S1H2 — 
molecule does not significantly change the HLG (6.9Sev) but on the other hand replacing Si with Sn significantly 
reduces the HLG (6.61ev). This is probably due to the metallic nature of Sn. 

We also carried out DFT simulations on the TP — SiH2~ molecule with Cu — S metal atoms replacing the terminal 
H atoms. The results indicate that the incorporation of the Cu atoms decreases the HLG energy gap to 6ASeV due 
to the transfer of charge from Cu atoms to the molecule. 

CONCLUSIONS 

Motivated by recent interest in the development of alternative to traditional electron devices, we have investigated 
a oligo-phenyl based double barrier molecular device. The double barriers have been created by using insulator 
spacers —CH2-, —SiH2— , GeH2 and SnH2 inserted between two aromatic rings We have carried out analytical 
as well as numerical studies using HF and DFT methods. We found that when the insulator spacers are same the 
tunneling matrix elements of the barriers are identical and this leads identical current flowing in both directions of 
the molecule. The symmetry can be broken by making the two insulator spacers different. This leads to asymmetrical 
electron conduction along either direction of the molecule leading to the desired characteristics of a typical diode. 
The asymmetric conduction properties such as the one investigated in this study are relevant for the implementation 
of molecular tunneling diodes. 
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